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Abstract 

A scheme which allows to compute the dynamics of strongly correlated classical ions embedded 
into a partiallzy ionized quantum plasma by first principles molecular dynamics is presented. The 
dynamically screened dust approach of Joyce and Lampe [Phys. Rev. Lett. 88, 095006 (2002] 
ist generalized to quantum systems. The electrons are treated fully quantum-mechanically taking 
into account their dynamical screening of the ion-ion interaction in linear response on the basis of 
an extended Mermin formula. The scheme allows to include the effect of the electron dynamics, 
electron streaming, wake effects and electron magnetization. 

PACS numbers: 52.65.Yy, 52.27.Gr, 05.30.-d, 05.60.-k 
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I. INTRODUCTION 



Strong correlation effects in ensembles of charged particles are of high importance in 
many fields of physics, including plasmas, the electron gas in solids or electron-hole plamas, 
and references therein. In recent years charged particles spatially confined in 



trapping potentj^i! 
dusty plasmas 



s have attracted considerable interest. Examples are ultracold ions 3, 41, 

1 n 

6|], or electrons in quantum dots [7], see Ref. m for an overview. 



II. BASIC EQUATIONS 

The system hamiltonian is given by 

H = J2Ha + J2Hab, (1) 

a ab 

= + W:\ (2) 



(3) 



HT' = E (4) 
Hab= E K6(|r.-r,|), (5) 

l<i<Na l<j<Nt 

where A and denote the vector and scalar potential of an external electromagnetic field. 
We will consider situations where the heavy particles, i.e. ions and neutrals are classical. 
Extensions to weak ion degeneracy will be briefly discusseed below in Sec. IV E[ Also, the 
neutrals will be assumed to have a minor effect on the ions which can be neglected at this 
stage. The effect of the neutrals on the electron dynamics will be included by an effective 
collision frequency, Uen = J2k ^en^ summing over all species of neutrals. 

The nonequilibrium dynamics of electrons (e) and ions (i) can be described by reduced 



density operators and is given by 



17| 



'^^t - [Hail), Fail)] = E^6Tr2[K,(l,2),F,,(l,2)], (6) 

where a = e,i, and Fabi^, 2) is the binary density operator and the operators are normalized 
to the volume according to |l7] 

Tri,..,Fi,.., = (7) 



. In the following, quantum exchange effects will be irrelevant and are neglected. Then we 
may introduce the binary correlation operator Qati^, 2) by separating the ideal contribution, 

While the ions are classical, no restrictions on the strength of their coupling are being 
imposed. In fact, we are particularly interested in strong ion- ion correlations leading to 
fluid and solid-like many-particle behavior. Therefore, the ion-ion interaction will be treated 
exactly. On the other hand, electron-ion coupling can be assumed weak, due to the quantum 
degeneracy (delocalization) of the electrons. Therefore, it will be sufficient to treat e-i 
interactions on the mean-field level, so e-i correlations will be neglected, Qei ^ 0. Then the 
ion equation becomes 
dp (1) 

'^—^ - W)^ = n^Tr2[Mh 2), ^,,(1, 2)] + n,Tr2[VUl, 2), F,{2)F,{1)] + 

n„Tr2[F,„(l,2),F,„(l,2)]. (8) 

The second term on the right gives rise to an effective mean field (Hartree potential) for the 
ions which is created by the electrons, Wie{l) = neTr2Vie{l, 2)Fe{2). The last term describes 
the interaction of ions with neutrals. 

Since the ions are classical we now take the semiclassical limit in the equation ([8]), first 



transforming to the Wigner representatio n 
by the Klimontovich phase space density 



ITj . Then the ion density operator is replaced 



N 



F, ^ Ni{r, p,t) = Y, ^[r - rfe(t)]5[p - Pk{t)] + 0{h^), (9) 



k=l 



which still includes all fluctuations. This function obeys a Liouville type equation 

{I + £s - ^ + + + ^-^^i 1} p- = °- ('°' 

Here (j)^^^ is the electrostatic potential due to all external fields and Win and Wie are the 
potential energies induced by the neutrals and electrons, respectively. Wie is obtained after 
an ensemble average of Wie over the electron density operator, i.e. Wie = {Wie)\F^ = Cj^e- 
The potential 0e is the solution of Poisson's equation and determined by the electron density 
via 

Mr,t) = ejd\'^^. (11) 

Furthermore, denotes the total potential energy due to all pair interactions of the 
ions. The superscript "Q" indicates that ionic quantum effects (which appear only at very 



short pair distances) have been taken into account by replacing the bare pair potential by an 
effective quantum pair potential, V^'^ . In the case of Coulomb interaction is the improved 
Kelbg potential, see Sj and references therein. Here we will, instead, have a pair interaction 
which is dynamically screened by the electrons giving rise to a different quantum potential, 
see Sec. IV E[ A similar procedure may be applied to the external potentials in the case that 
they vary on small length scales comparable to the size (or thermal DeBroglie wave length) 
of an ion, e.g. [20^, but we will not write this explicitly. 

Eq. (fTOj) is equivalent to a set of Newton's equations for the ions (which are its charac- 
teristics) which is directly verified by introducing the definition ([9]) into (fTOj) 

^ = -^(^^ + + e^<^"*)|r=r, - ''^nPk + y{t), k = 1, . . . , N„ (12) 

which treat the ion interactions and the external fields exactly. In case of a time- dependent 
electromagnetic field, the r.h.s. will contain the ionic Lorentz force instead of the gradient 
of 0*^^*. The last two terms on the r.h.s. are related to the potential Win in Eq- ( |TOll and 
allow to simulate the effect of the neutrals on the ions (friction term and Langevin noise 
term), but may be neglected in most cases, except for a verly low degree of ionization. In 
addition, the effect of the electrons is taken into account via the induced potential (pei^jt), 
Eq. (ITT]) . The quality of this potential depends on the level of treatment of the electronic 
density which is computed from the electron density operator. We, therefore, now turn to 
the computation of F^. 



III. QUANTUM ELECTRON DISTRIBUTION FUNCTION IN LINEAR RE- 
SPONSE 

We now return to the equation of motion ([6]) for the electron density operator. Since 
we consider a dense degenerate electron system with small value of the electron Brueckner 
parameter, r^e = r^/aB < 1, electron-electron correlations are small. The main contribution 
to the interaction between electrons will, therefore, arise from the electronic mean field. Wee, 
which will be treated exactly. The additional correlation effects between electrons and from 
scattering with neutral particles (contained in the electron binary correlation operator) will 
be treated approximately, within a relaxation time approximation. This term is denoted as 
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collision integral /e, its form will be specified below. 



ih^^-mi),F,{l)] = /,(!), (13) 
H^il) = Hl{l) + #ee(l), (14) 

W,e{l) = neTr2Veeil,2)Fe{2). (15) 



A. Momentum representation for a spatially inhomogeneous system 

In the following we will consider the response of the electron system to a spatially inho- 
mogenous electrostatic field. To do this, it is convenient to transfrom the abstract operator 
equation (fT5|) into the momentum representation using as basis function free electron states 
\k). Thus, we split the single-particle hamiltonian into a field-free part, h^, and a field part, 
U, = hl + U, with 

hl\k) = ek\k), ^it\iek = {hkf/2me. (17) 

In this basis, the operators transform into matrices (we drop the time-dependence): 
{k\Fe\k') = V~^fk,k', (^l^el^') = ^kSk,k', {k\U\k') = Uk,k'- The normalization condition 
2 fkk{t) = Nf,{t) where the factor 2 arises from the spin summation. Multiplying Eq. 
( ITSl) from left by {k\ and right by \k') it becomes a matrix equation 

^h^^ - (efc - ek')fk,k' - (fc|[f/^^, FJIA;') = h,k', (18) 

where the effective potential operator is U'^^ = U + W^e- The evaluation of the commutator 



term will be performed below in linear response. For the full result, see Ref. 17|. 



B. Electron equation in linear response 

We will now linearize Eq. (fT8|) for the case of a weak external potential JJ around the 
field- free solution /o, 

/fc.fc' = fk,k' + ^fk,k', with \5fk,k\ « fk- (19) 
Without external field it obeys a spatially homogeneous equation, i.e. /°^, = fkSk,k', 

= II (20) 
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where the other terms in Eq. (fTS!) vanish due to homogeneity. In the colhsionless case, 
the unperturbed solution is taken as the equihbrium distribution - in the quantum 
Fermi distribution, 

/f^=[e^(--^) + l]-\ (21) 

where /3 = l/ksT and fi{n,T) is the chemical potential. This is correct for Markovian 
collision integrals which vanish for Fermi-Dirac functions, i.e. Ikifk'^] = 0' t)ut neglect 
nonideality effects 17|. In principle, instead of Eq. (12T|) any other stationary solution is 
possible. For generalized non-Markovian collision integrals one first has to solve for the cor- 
related distribution function, as it can be done e.g. within nonequilibrium Green's functions 
theory, e.g. [24 1. 

The next step is the solution of the linearized equation for Sfk,k which is driven by an 
external electromagnetic perturbation and, in the coUisionles case, gives rise to the Lindhard 
(RPA) dielectric function (DF) or, in the classical limit, to the Vlasov DF. In case of a 
monochromatic excitation the perturbation of the distribution oscillates resonantly with the 
field without irreversible relaxation. There is no back reaction on the function fk- 

Here, we want to go beyond this approximation and include the effect of collisions into 
the dielectric properties and in the induced electrostatic potential. Thereby we will assume 
a Markovian collision term which allows us to use Eq. fl2Tl) for the unperturbed solution. 
We first write down the equation for the perturbation 6fk,k which is obtained by linearizing 
Eq. ( fT9l) around the homogeneous field-free solution fk, 

^^ ^^/fc.fc _ _ _ ^f/cff^ . (j^, _ ^ 6Ik,k', (22) 

lim fk,k'it) = fkSk,k'- 

t—t — OO 

Here, 6U^^ = U + 6Wee is obtained by replacing f^^k' by 6fk^k', whereas 61 is obtained by 
kee pin g in all appearances of the electron distribution functions only terms of first order in 
Sf 2J]. Due to linearity, this equation is conveniently analyzed in Fourier space. Assuming 
a monochromatic perturbation, U{r,t) = Uq{uj)e^^^'^^^^^^~^^'^^ the same time and space de- 
pendences occure in all terms. Here S is an infinitesimal positive real number assuring that 
the potential vanishes at t ^ — oo. 

To relate the wave vector q of the spatial modulation to the momenta k, k' we introduce 
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center and relative momenta by 

Q = — - — , q = k — k , or, vice versa, (23) 
k = Q + l, k' = Q-^-. 

Then, Eq. ( !22l) can be written in Fourier representation as 

[h{uo + i6) - (efc - e,0] ^fk,k' - SUf{uj){h' - h) = SI,{uj), (24) 

where the common exponential factor has been cancelled. 

C. Collision integral in relaxation time approximation 

If collisions are present they give rise to relaxation to a stationary state which will be 
achieved for times exceeding the relaxation time r of the system. In the presence of a 
monochromatic excitation, this asymptotic solution will be modulated by the field and will 
not conoincide with the Fermi function (pT]) . It can be constructed as a sum of the stationary 
solutions of Equations fl20l) - i.e. the Fermi function - and the asymptotic solution of Eq. 
(l22|l . The latter is readily obtained by putting the time derivative to zero and requiring the 
collision integral to vanish, SIk^k'[^f°°] = 0. Correspondingly, in the Fourier representation 
we consider the limiting case u —* and obtain the stationary asymptotic solution: 

Sf^,, = ^^^SU^,r, co = 0. (25) 

The superscript "oo" of the effective potential indicates that the asymptotic solution 6f°° 
has to be inserted but also that the static limit of the effective potential has to satisfy 
an additional consistency condition. As we will see below in Sec. IIIIDI this is related to 
fulfillment of sum rules and conservation laws. 

Thus we found the asymptotic solution of the full electron distribution function 

lim fk,k'(t) = fkSk,k' + Sf^k'^^fl (26) 

and, for 7^ 0, 

lim fk,k'it) = fk6k,k' + 5fkA^y^'- (27) 

t— *+oo 

Since 5f is a small correction, the stationary result (!26|) can, in linear order, be 
included into the Fermi function giving rise to a local Fermi distribution f^^{r) = 

7 



^Q0iik-fj.-Sfi{r)) _|_ 1^ where Sfi{r) = 6U°^'°°{t) [21], consistent with thermodynamic sta- 
bihty in an external field, but this is not necessary for the derivations below. 

After having found the asymptotic solution 5/^/ of the distribution function we can 
construct the collision integral 61 which drives the relaxation towards this result. A simple 
static model which has this property is the standard relaxation time approximation. We 
recall the result for the spatially homogeneous case (we drop all arguments), 

| = /^^^--^(/-r'5), /(0) = /o. (28) 

The solution of this equation together with the initial condition at t = is f{t) = /oe^*/"^ + 
f^^[l — e^*/"^], showing the decay of the initial state and the approach to the asymptotic 
state. Using this result we can now construct the collision integral 61 of Eq. fl2^ in RTA. 
Replacing in Eq. fl28|) / by 6f and f^'^ by 6fj^j., [the function cancels] and multiplying 
by we obtain 

T I Cfc - Cfe/ 'J 

Here, r is the relaxation time which has to be computed from a separate kinetic theory 
or taken from experiment. We use a simple static approximation where r is frequency 
independent which is consistent with the asymptotic result (|25l) which was obtained for 
u = 0. Recall that 6Ul^f^^ still has to be determined. 



D. Distribution function 6f with collisions 

Using the result ( 129|) in Eq. ( 12^ we are now ready to explicitly compute the perturbation 
of the distribution function. A straightforward calculation gives 

{i-TTcfTjOO "\ re 
6Ul% {u) - rhv^^ J.\-\' , (30) 

where we used 5^0 due to the existence of a finite collisional damping v = t~^. This 
result is a straightforward extension of the collisionless random phase approximation. Scat- 
tering effects (terms proportional to z/) are contained in two places: first, the frequency in 
the denominator is replaced by a complex frequency and, second, there appears an addi- 
tional contribution proportional to 6U^^'°° in the numerator which renormalizes the Fourier 
component of the effective potential. We now turn to the computation of this latter term. 
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To this end we consider the local particle conservation law, i.e. the continuity equation 

^^ + divj(r,t)=0, (31) 

where density and current density of the electrons are related to the distribution functions 

by 

nir,t) = 2 J ^^/,(p,r,t), (32) 

and the prefactor 2 accounts for the spin summation. Now, in the field-free case, the 
distribution function is time- independent and spatially homogeneous and does not contribute 
to the continuity equation. In contrast, in the presence of the external field U there is a time 
and space dependent contribution 6f which contributes to n and j. Thus, the continuity 
equation becomes, after transformation to Fourier space, 

uj6nq{uj) = q ■ 5j<j(cj), (34) 

where 

5ng{uj) = 2 J ^^-^5/q+|,q_|(^), (35) 

W = 2/||^.fc,,«_,(.). (36) 



We will now verify that local particle conservation (IMI) is satisfied by our solution which 
gives the required condition for if^^'^". To compute the density and current from the solution 
fl30|) requires to perform integrations over Q. To shorten the notation we introduce the 
following relevant integrals (n = 0, 1, . . . ) 

n^,n = 2 / j—^ — X — — -, (38) 

where we introduced the abbreviation for the complex frequency Co = uj + iu. In the colli- 



sion 



ess limit, u —>■ 6, and the n„ become (moments of) the retarded polarization function 
Further, we introduced the short notation Q±^ ±. Note that the integrals with 
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odd powers are vectors in the direction of Q. We these definitions the Fourier components 
of the density and current density become (due to homogeneity, SU^j., = SU^^) 

Sn,{uj) = Uo{q,u;)SUf{uj)-thiySUf'°°U,,o{q,uj), (39) 
6j,{uj) = U,{q,u;)6Uf{uj)-thiy6Uf'^U,,i{q,uj). (40) 

We now transform times the integral Hi, by adding and subtracting under the integral 
h(2}. Taking into account that e_|_ — e_ = h^Q ■ q/mg, we obtain the identity 

hciIl,{q,Lu) = -2 J ^(/- - /+) + {nLU + thiy)Uo{q,u). (41) 

Assuming that the field-free distribution depends only on the modulus of the momentum, 
i-e. /-k = /k the integrals over /_ and cancel. The same transformation is possible for 
the integral with the result 

hqUi,^i{q, u) = Uoiq, 0) + {hu + ihu)U^^o{q, i^)- (42) 

Collecting the results (HTj) and P2|) together we may rewrite the expression for the current 
density, Eq. (HOj) 

hqSi,{uj) = h{uj + iiy)Uo{q,u:)SUf{uj) 

-iTiv [no(g, 0) + [Tiuj + i;iz/)n,,o(g, u:)] 6Uf'°° = 
= hudnqiuj) + 

+ihu {no(g, uj)6Uf{uj) - t/f [no(g, 0) + 2^J/n,,o(g, uj)] } . (43) 

Evidently, the continuity equation (!34l) is fulfilled if the terms on the last line (in the curley 
brackets) vanish. From this we find the condition for the asymptotic value of the effective 
potential 



q TT ^„ I „-fe,.TT f„ ,^:\ TT f„ r\\' \ ' 



Yio{q,Qj)5Uf{uj) _ 5n,{u) 
no(g,0) + i^z/n,,o(g,cI') no(g,o)' 
To obtain the last equality we used YlQ^q, Cj)5U'^^ {u) = dnqioj) + ihuYiy^Q{q,Cj)U^^'°^ . With 
this the problem has been solved and the perturbation of the electron distribution function 
(p7]) has been obtained explicitly: 

TeSf..^ ^nq{uj) \ fk' - fk 



no(g, 0) [efc - ek'\ J n{uj + iv) - [e^ - ty] 
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IV. QUANTUM DIELECTRIC FUNCTION CONTAINING COLLISIONS 



To compute the dielectric function we first have to obtain an exphcit result for the density 
fluctuation 5n. Using fl5I?l) and (jUj) or integrating fH5]) . we can write 

5n,{u) = Uoiq,co)SUfico)-thu^^^^^Sn,ico). (46) 

Recalling the definition of the effective potential, 6U^^{uj) = 5Uq{uj) + Vq5nq{(jj), and intro- 
ducing the short notation Yiy Q^q^Cj) = Iiy^Q{q,Cj) /liQ^q, 0) we can solve for 5nq\ 

5nq{uj) = ^^il^^) ^ ^jj^^^y (47) 

1 - Vgno(g, uj) + ihuU^fl{q, to) 

The inverse dielectric function is defined as |l 

and is obtained by inserting the solution fl47l) into the effective potential with the result 

6Ufiu) = l + ^^^n^.o(^.^) (49) 

1 - VgUo{q, uj) + ihuIl^^o{q, u) 

Performing the derivative with respect to 6U and inverting the result, we obtain the dielectric 
function 

1 + ihi'Ll^^o[q,uJj 

An alternative way of writing this result is to eliminate the function IIi^^o by using the identity 
hujll^fi{q,u}) = Ilo{q,uj) — no(g,0). In this case, in the expression for the dielectric function 
one can make the replacement Il^^o{q,uj) = [no(g, u))/no(g, 0) — l]/hu!. 

This way we have succeeded to derive a dielectric function which contains collision effects 
in a relaxation time approximation with a static collision frequency u. This is the result of 
Mermin j2l| which had been derived before for classical plasmas by Rostoker and Rosenbluth, 



e.g. 







A. Further improved dielectric functions 

Further improvements of the Mermin result have been considered by various groups. 



Ropke et al. have derived a 
contains energy conservation 



VIermin-type expression which, besides particle conservation 



27|. However, they found that the effect was small. Another 
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modification by this group was to include a frequency dependent collision frequency into 



the relaxation time collision integral 



Finally, we mention that a selfconsistent nonequi- 



librium calculation within Nonequilibrrium Green's functions which fully included sum rule 



preservation has been recently performed 



B. Dynamical screening and wake effects 

The main motivation to include dynmical screening of the interaction between heavy 
particles in a two-component plasma is its importance for nonequilbrium situations. One 
such case is the existing of streaming light particles which causes wake effects which have 
a dramatic effect on the arrangement of heavy particles, e.g. dust particles in a complex 
plasma. This was discussed in detail by Joyce and Lampe, cf. 22] and references therein. 
Wake effects in a quantum plasma have also been considered by one group [25|| who found 
an important influence on stopping power of ions in a polarizable medium. However, these 
were only single particle effects. 

V. DYNAMICALLY SCREENED ION-ION PAIR POTENTIAL 

Let us now compute the potential of a moving charged particle taking into account the 
dielectric properties of the plasma. We start with the case of a classical charged particle, 
e.g. and then generalize the result to quantum particles. 

A. Potential of a moving classical particle 

The Poisson equation for a polarizable medium reads 

divD(r, t) = CaUair, t), (51) 

which in Fourier space reads ih ■ Dfc(t) = eank{t). The electrostatic potential created 
by the charge density on the right is E(r) = — V0(r), corresponding, in Fourier space 
to Efc = — ik ■ (pk- Together with the electrodynamic definition of the dielectric tensor, 
Dk,i = Zlj ^k,ijEk,j, i,j = 1, 2, 3 and Eq. (EI]) we obtain 
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Consider now the case of a classical point charge "a" with initial position r(0) = roa, moving 
with constant velocity Vq (relative to the carriers creating the dielectric function). Then 

na{r, t) = 5[v - - Vat], (53) 

with the Fourier representation 

nakiio) = 2ne'''"^6[iu-ky^a]. (54) 
Inserting this result into ( l52l) and performing the back transform we obtain 

/d'^k g«k(roa-r) 
/ / n ^• (^^) 
i^^r 22ij k^kjek,^j{k^r^) 

Thus, the potential of a classical point particle moving with a constant velocity is time- 
independent, only a single frequency component, u = kv^, is present in the spectrum. 
For the special case of an isotropic medium, e^j has only two independent components. 
Concentrating on longitudinal plasma oscillations we can replace eij kikje/k"^, and the 
potential (!55|) becomes 

/d^k e e*'''-'"'''""'"^ 

This result may be immediately generalized to the case of many particles. Indeed, due 
to linearity of Maxwell's equations, the resulting total potential is simply the sum of all 
potentials of the type i.e. 0*°*(r) = Y.a=i ~ ^Oa] v^). 

B. Potential of moving quantum particles 

In case the charge density is created by quantum particles which are not point like or by 
an ensemble of many particles the particle density is expressed, e.g. via the single-particle 
distribution function. 

On the other hand we can consider the result (!55l) in the continuum limit of a superposition 
of infinitesimal charges. Then the contribution d(f) from the charge dq{r) inside a small 
volume dV around point r is proportional to dq{r)/dV = en{r), and the total electrostatic 
potential is 

J (27r)2 F efc(kv,)' ^ ^ 
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where n has to be computed from the distribution function according to (157|) . and it has 
been assumed that all particles stream with the same velocity v^. This result is applicable 
to classical and quantum systems equally and provides, in particular, the induced potential 
of an ensemble of degenerate electrons with e characterized by a distribution function f^.. 

C. Ion-ion pair potential 

We can now return to the ion djTiamics studied in section [III There the effect of the 
electrons was comprised in the electrostatic potential 0e, cf. Eq. f lTUl) . If the electron density 
is spatially homogeneous, the particles will produce a homogeneous potential 0(r) = const, 
as can be seen from Eq. (!58|) . and will not exert a force on this ions. In an external potential 
the electron density maybe spatially modulated but this will be a small effect. Therefore, the 
main effect of the electrons is polarization of the plasma medium via the dielectric function 
( l50l) . In this case, the electrons will renormalize (screen) any other potential created by 
charged particles, in particular the ion-ion interaction, term in Eq. (ITOl) will be screened, 
i.e. — s> . Due to the linear approximation, we may write 

V?{v)= (59) 

l<k<l<Ni 

The corresponding force term in the Newton's equation for the /c-th ion will be minus the 
gradient of 

ek4>tot{rk) = ^ (rfc - r;) = Cfc ^ (j)i{rk-ri). (60) 

l<l<Ni,l=/=k l<l<Ni,l=/=k 

The screened potential created by the l-th ion which moves with velocity Vi(t) (relative to 
the electrons) then follows directly from our result fl56|) . Note the direction dependence of 
the potential and the dynamically screened pair potential. If the electrons stream with a 
constant velocity Ug this potential becomes 

/(i^k 1 e«k(roi-r) 
(27r)^ Cfc (k ■ [V;(t) - UeJ) 

and may be time-dependent through the time dependence of the velocity vi. Further it is 
straightforward to include ionization and recombination effects through a time dependence 
of the charges Cfc and e;. 
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D. Limiting cases 



The classical limit of the potential (16T|) gives the Rostoker-Rosenbluth result 23| which 



was used in classical MD simulations of dusty plasmas [22|. The present potential is a 
straightforward generalization to dynamical screening by quantum particles. The relevance 
of dynamical screening depends on the particle velocities. If the electrons are at rest, we will 
have, on average Ug = 0. If further, the ions are at rest or moving very slowly, the frequency 
argument of the dielectric function vanishes and k'^/{k'^ + k^). In thermodynamic 

equilibrium the screening parameter is given by the Thomas-Fermi length, ^ vtf- Then 
the potential (1611) contains just the Fourier transform of + which yields the familiar 
Yukawa potential. Thus the resulting ion-ion pair potential reduces to the isotropic static 
Yukawa pair potential Vii{r) e'^e~'^^/r. Thus the validity of the Yukawa approximation 
is limited to electrons at rest and to the neglegibility of ion thermal motion as well as of 
quantum degeneracy effects. 



E. Quantum effects in the ion-ion pair potential 

Let us briefly comment on the quantum effects in the pair interaction (denoted by the 
superscript "Q"). In using the expression fl56|) we assumed pointlike ions and neglected any 
finite ion extension. This can be corrected in the final expression by replacing the Fourier 
transform of the Coulomb potential, i.e. the factor l//c^ in Eq. fl6Tl) . by the Fourier transform 
of the improved Kelbg potential [8] or of any other appropriate quantum potential. 



VI. DISCUSSION 



A new model for the simulation of dense quantum plasmas including dynamical screening 
of the electrons, partial ionization and strong ion correlations has been developed. It is par- 
ticularly important for high density low-temperature plasmas, in situations where the ions 
form liquid or solid-like structures and when the electrons are in nonequilibrium. Typical 
situations are streaming electrons or plasma instabilities due to fast electrons or electro- 
magnetic fields. The proposed simulation scheme is based on classical molecular dynamics 
simulations where the dynamical screening effects are incorporated within a linear response 
approach into the screening of the ion-ion pair interactions. For the screening an extension 
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beyond the (collisionless) RPA model has been used which is due to Mermin and a strict 
derivation within quantum kinetic theory has been given putting the original result of Ref. 
2l| on solid ground and critically assessing its scope of applicability. 
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